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Abstract 

Folded supersymmetry (/-SUSY) stabilizes the weak scale against radiative cor¬ 
rections from the top sector via scalar partners whose gauge quantum numbers 
differ from their Standard Model counterparts. This non-trivial pairing of states can 
be realized in extra-dimensional theories with appropriate supersymmetry-breaking 
boundary conditions. We present a class of calculable /-SUSY models that are 
parametrized by a non-trivial twist in 5D boundary conditions and can accommodate 
the observed Higgs mass and couplings. Although the distinctive phenomenology 
associated with the novel folded states should provide strong evidence for this 
mechanism, the most stringent constraints are currently placed by conventional 
supersymmetry searches. These models remain minimally fine-tuned in light of LHC8 
data and provide a range of both standard and exotic signatures accessible at LHC13. 



2 


I. Introduction 

The canonical solutions to the naturalness problem of the Standard Model are based on 
supersymmetry (SUSY) or compositeness. Both of these approaches lead to “top partner” 
states that carry the Standard Model gauge quantum numbers of the top quark and couple 
to the Higgs boson with a strength set by the top Yukawa. The QCD charge of the top 
partners implies a large production cross section at the LHC. Searches for these states have 
been performed using data from LHC8, and the null results can be interpreted as constraining 
naturalness in the context of these models (see e.g. [1, 2] for reviews on the status of SUSY). 
In particular, the non-observation of partner states gives rise to a “little” hierarchy problem 
- percent-level cancellations between threshold corrections are required to reconcile the 
weak scale with experimental bounds on the new states predicted by supersymmetry or 
compositeness. 

However, as was first pointed out in the seminal papers [3, 4], it is logically consistent 
for the top partner states to be color-neutral. The resulting signatures can differ wildly 
from those traditionally associated with naturalness. The production cross-section of the 
un-colored top partners is significantly reduced, with correspondingly weaker limits at the 
LHC [5] . Lighter top partners are thus consistent with data, ameliorating the tension of the 
little hierarchy problem. There are two broad classes of these so-called “neutral natural” 
models: Twin Higgs models [3, 6-11] (or their generalization into the Orbifold Higgs [12]) 
rely on a shift symmetry to protect the Higgs potential, whereas models of Folded SUSY 
(/-SUSY) are supersymmetric at short distances [4], In both cases, discrete symmetries 
enforce the equality of the Higgs couplings to the top quark and its un-colored partners. 
This results in the accidental partial restoration of the UV symmetry at low energies. 

For models with shift symmetries in the UV, e.g. Twin Higgs, the Higgs quartic must 
necessarily be generated at loop level (with the possible exception of terms arising from 
collective symmetry breaking [6, 13]). Since the Higgs mass parameter is generically one loop 
down from the scale of global symmetry breaking /, naturalness suggests that v ~ /, where v 
is the Higgs vacuum expectation value (vev). Furthermore, misalignment of the Higgs boson 
from the SM vacuum leads to irreducible Higgs coupling deviations parametrized by (v//) 2 ; 
experimental constraints derived from Higgs properties require / > 700 GeV [14], imposing 
a minimum tuning cost irrespective of the direct search reach for partner states. In contrast, 
neutral natural models relying on supersymmetry, such as /-SUSY, can accommodate both a 
tree-level quartic and a Higgs boson aligned with the Standard Model vacuum. As such, top 
partner masses can be pushed as high as ~ 500 GeV with no attendant loss of naturalness. 

How does one construct a model of neutral naturalness that is consistent with SUSY 
at high energies? Specifically, since the J\T = 1 global supersymmetry commutes with the 
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Standard Model gauge group, the low-energy pairing of colored and un-colored states in /- 
SUSY resembles hard SUSY breaking from a 4D perspective. Clearly both the colored 
and un-colored states must come in full supermultiplets. Then the desired low energy 
spectrum must be realized through SUSY-breaking. This can be naturally realized in 5D 
theories compactified on S i/Z 2 with suitable boundary conditions. Since SUSY breaking 
is realized non-locally in the extra dimension, i.e., the Scherk-Schwarz mechanism [15, 16], 
contributions to the Higgs potential are finite despite the inherently low cutoff scale typical 
of 5D models. 

Given this feature of calculability, it is interesting to undertake a careful analysis of 
electroweak symmetry breaking in /-SUSY. In this paper we describe the relation between 
the allowed parameter space and weak scale fine-tuning. We will begin by analyzing /- 
SUSY in its original incarnation, wherein the un-colored /-top partners are massless at 
tree-level while their colored counterparts have a large mass set by the compactihcation 
scale. As will be shown below, the Higgs potential (given by an approximate two-loop 
calculation) does not appear to yield a non-zero vacuum expectation value for the Higgs 
boson. This motivates deforming minimal /-SUSY in order to simultaneously achieve both 
electroweak symmetry breaking and a 125 GeV Standard Model-like Higgs boson. One 
minimal extension is to allow for a non-trivial twist a G [0,1/2] in the Scherk-Schwarz 
boundary conditions, leading to tree-level masses for the zero-mode /-scalars and thereby 
controlling a partial cancellation between top and /-top sector contributions to the Higgs 
potential. This additional parametric freedom provided by a allows electroweak symmetry 
breaking to be triggered at one loop without spoiling the attractive features of finiteness 
and calculability. 

In the minimally tuned regions of parameter space, /-SUSY models with non-trivial 
twist give rise to a range of observable phenomenology. Generic predictions include many 
of the exotic signatures associated with the original model of /-SUSY. For example, the 
/-squarks with mass \a — l/2|/i? behave as scalar squirks [17] that annihilate to electroweak 
resonances [5,18]. There are also interesting hidden-valley-like [19, 20] observables associated 
with the /-glueball states [11, 21, 22], Additionally, these models may be probed indirectly 
through precision Higgs measurements [23, 24]. 

The precise phenomenology depends in detail on both the Scherk-Schwarz twist and 
on the geographic distribution of supersymmetric multiplets in the extra dimension. For 
concreteness, we will focus on two specific models distinguished by the details of the Higgs 
sector. In the first model, the Higgs bosons will be restricted to a four-dimensional brane. 
In this case, we will see that the dominant constraint arises from a combination of precision 
Higgs measurements and searches for the heavy Higgs bosons. In the allowed region of 
parameter space, this model is tuned at the 0(0.2%) level. This motivates studying a second 
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model where the Higgs fields are allowed to propagate in the bulk. Now both the Higgsinos 
and the heavy Higgs states arc lifted using boundary conditions. This alleviates the tension 
with Higgs measurements and allows for a reduction in fine-tuning to better than the 10% 
level. The dominant collider signature will come from gluino-squark associated production. 
Intriguingly, compressed spectra are a generic feature of these models, weakening limits on 
radius of the extra dimension. In both the brane and bulk Higgs cases, discovery will proceed 
through the detection of MSSM-like states, followed by the spectacular observation of the 
/-squark sector. 

The rest of this paper is organized as follows: In Sec. II we review the generic features of 
5D supersymmetric theories with Scherk-Schwarz SUSY breaking. In Sec. Ill we describe 
the top and gauge sectors of /-SUSY with arbitrary Scherk-Schwarz twist, along with 
the radiatively-generated Higgs potential. In Secs. IV and V we present in detail the 
parameter space and fine-tuning implications of two complete models in which the Higgses 
are alternately localized on the brane at y = 0 or allowed to propagate in the bulk. Finally, 
we conclude in Sec. VI. We reserve a variety of technical results for a series of appendices. 

II. Supersymmetry in 5D 

In this section we review some basic properties of theories with Scherk-Schwarz SUSY 
breaking. In particular, we consider 5D theories compactified on S\ /Z 2 in which SUSY is 
broken by boundary conditions on the orbifold fixed points. As the smallest spinor in 5D 
is Dirac (with 8 real components), when written in terms of 4D superfields [25] the matter 
content must be equivalent to that of = 2 SUSY. Furthermore, the interactions in the 
bulk must respect the full SU(2) x U (1)r symmetry. Just as in the non-supersymmetric case, 
the low-energy theory will not contain chiral fermions unless suitable boundary conditions 
are chosen to “project out” the zero-mode Dirac partners. A Dirichlet (Neumann) boundary 
condition then corresponds to an odd (even) eigenvalue under reflection about y = 0 (tt R), 
hereafter denoted Z (. Z '). This may be accomplished in an A^ = 1 SUSY-invariant way by 
choosing Z and Z' to act on complete 4D superfields. 

A. Supersymmetry Breaking by Boundary Conditions 

The action of the transformations Z, Z’ on superfields are 

Z[*{y)\ = Z${-y) Z'[<t>{y-*R)\ = Z' ${-y + tt R), (1) 
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where Z 2 = Z /2 = 1. The apparent M = 2 SUSY is broken down to A/" = 1 if the eigenvalues 
are different for each of the Af — 1 superfields in an A^ = 2 multiplet. In particular, we 
may choose Z' = US 3 [/ _1 , in which S 3 is +1 (—1) on even (odd) M = \ superfields. Here 
we take U to be a 2 x 2 matrix acting on the SU(2)r doublets; it acts trivially on SU(2)r 
singlets. More generally, the eigenvalues E 3 can be chosen to break additional symmetries, 
and likewise the rotation matrices can act on other symmetry multiplets [26], although we 
ignore this possibility for simplicity. 

We will work in the basis in which Z is diagonal, Z = S 3 . If boundary conditions are 
chosen such that the same Af = 1 SUSY is respected by both of the fixed points, then 
U = 1 and Z' = S 3 , corresponding to [Z, Z'\ = 0. In this case, the low-energy theory 
after KK decomposition will also be AT = 1 supersymmetric. In general, however, the fixed 
points may be taken to be invariant under inequivalent SUSY transformations. Then we can 
parameterize this commutator by a single angle (ignoring possible phases), 0 < 2na < n, 
which parameterizes the rotation between bases in which Z,Z r are respectively diagonal. 1 
In the diagonal basis of Z, Z' is then given by 


_ g— ■Kio.ci ^7 g7ria<72 


( 2 ) 


For a = 0, Z and Z' commute as above, and the theory is AS = 1 supersymmetric. At 
a = 1/2, these Z 2 transformations are “orthogonal” in the sense that again [Z, Z'~\ = 0, but 
Z' = —Z. For general 0 < a < 1/2, Z and Z' do not commute; in the basis where Z = S 3 , 
we have 


[Z,Z'\ 


^2 7T i ex. cr 2 g —2 ni ot <72 


( 3 ) 


As an example, consider a matter hypernmltiplet, (<&, $ c ), composed of a conjugate pair 
of 4D chiral superfields, $ ~ (0, -0) and <F C ~ (</> c , -0 C ). Without loss of generality, we may 
take $ and <F C to be eigenstates of Z. Then the eigenstates of Z' are, in general, inequivalent 
chiral superfields < F / ~ (//, U) and <F C/ ~ (0 C/ ,'0 C )> in which 




( 4 ) 


At maximal twist (a = 1/2), </Y = /U and 0 C/ = —cjA. This situation is illustrated in Fig. 1. 
After integrating out the extra dimension, the component fields may be KK decomposed, 
ip(x, y ) oc Yhk fk(y) ^ k \x)-, with each mode having an apparent 4D mass m\ fk = dy fk- For 
example, given the boundary conditions in Fig. 1, the KK modes of each component field 


1 If other symmetries were also broken with boundary conditions, the commutator would be parameterized 
by additional angles corresponding to the orientation of Z , Z' in the space of these additional symmetries. 
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FIG. 1: Component field eigenvalues under reflections at y = 0 and y = it R, shown for a matter 
hypermultiplet (left) and a vector multiplet (right) at maximal twist. Solid (dashed) lines represent 
component held groupings into 4D superfields at y = 0 (ttR). 


will be split in mass. At low energies this is just a theory with softly-broken SUSY. This 
approach to SUSY breaking is the Scherk-Schwarz mechanism; for explicit wave functions 
and arbitrary boundary conditions, see the Appendix. 

Noting that the theory is locally J\T = 1 invariant everywhere in the extra dimension, all 
SUSY-breaking effects must be global in nature. In particular, all UV divergences, which 
correspond to contractible loops, must vanish supersymmetrically. The SUSY-breaking 
contributions of these loops come from wrapping the virtual propagators around the extra 
dimension, so that they “see” the SUSY-breaking boundary conditions. The space-like 
separation of these propagators gives an exponential suppression at high momenta, rendering 
all such loops finite [27-31]. Therefore, although 5D models with bulk gauge fields are 
intrinsically non-renormalizable and necessarily come equipped with a UV cutoff A (which 
in the models we will consider is not far above 1 /R), all SUSY-breaking effects will be 
calculable. In particular, as we will see in the following sections, all corrections to the Higgs 
potential are explicitly finite. 

Another generic feature of Scherk-Schwarz SUSY breaking is a compressed spectrum of 
all the 4D states who descend from bulk fields [32, 33]. All of the bulk states of the theory 
share a common A-symmetry; absent additional global symmetries which may be broken by 
boundary conditions, common boundary conditions for all of the superpartners are enforced. 
Thus, up to radiative and/or volume-suppressed corrections, there will be a common mass 
for all zero-mode superpartners. Of course, brane-localized states do not share this common 
mass, and so ultimately the implications of compression in Scherk-Schwarz theories depends 
on the detailed geography of brane and bulk states. As we will discuss in Section V, there 
are a variety of realistic models in which compression leads to a considerable reduction in 
the strength of collider limits. 
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B. Calculability and its Discontents 

The novel properties of Scherk-Schwarz SUSY breaking have been explored in previous 
papers in the context of the Minimal Supersymmetric Standard Model (MSSM) [26, 27, 
29, 33-46]. In addition to a natural, calculable breaking of electroweak symmetry, these 5D 
implementations of the MSSM predict a relatively large Higgs mass (quite near the measured 
value) when one or both of the top quark multiplets propagate in the bulk. However, since 
the scalar potential is controlled by only one scale, the compactification radius is uniquely 
predicted by the Higgs vev, and there is no parametric freedom that could allow the tuning 
of radiative contributions in order to raise 1/R while keeping the Higgs vev fixed. As 
might be expected from naturalness considerations, one then finds the observed weak scale 
corresponds to 1/R ~ 700 GeV [34], and therefore an experimentally-excluded gluino mass, 
rrig ~ 350 GeV. 

Introducing the folded sector provides the necessary parametric freedom to raise the KK 
mass scale beyond the reach of current bounds. In the limit of maximal twist (see Fig. 2 
below for an illustration) the one-loop top contribution to the Higgs potential vanishes due 
to the presence of an “accidental supersymmetry.” As a is lowered from a = 1/2, the 
size of the top contribution to the Higgs soft mass is increased from zero, while the gauge 
contribution is decreased, allowing the value of 1/R corresponding to v ~ 174 GeV to be 
continuously varied. In folded theories, the weak scale can be accommodated with larger 
values of 1/R relative to those predicted by 5D realizations of the MSSM. 

Alternatively, one could relax the requirement of calculability and couple the MSSM to 
arbitrary SUSY breaking on one of the branes; this possibility has recently been explored 
in [47]. In these models, there is a new scale of SUSY breaking that is not determined by 
the radius of compactification. In principle, electroweak symmetry breaking along with an 
arbitrary spectrum can be achieved, but at the expense of introducing irrelevant operators 
with relatively large, unknown coefficients. For the remainder of this paper, we will neglect 
this possibility by focusing on calculable scenarios in which the Scherk-Schwarz twist is the 
only source of SUSY breaking. 

There is one caveat to complete calculable control. Any operators allowed by AT = 1 
SUSY that are consistent with all the symmetries in the theory may be generated on the 
4D branes at the cutoff scale A. This manifests as non-trivial brane kinetic terms and other 
incalculable Kahler potential operators, such as K D 5(y) Z l0 dr,-. As the theory is strongly 
coupled at scale A, these coefficients should not be assumed to be small; in particular, they 
may have a substantial impact on the low-energy spectrum since A R is not large. 

One concern is that the flavor structure of these terms would lead to large flavor violating 
effects at low energies. Since this is set by the strong dynamics above A, we will assume 


that the UV theory respects the flavor symmetries of the Standard Model up to the Yukawa 
couplings, i.e., minimal flavor violation. We will therefore not consider flavor constraints 
on the model from large brane-localized operators. Of course, there may still be sizable 
flavor-preserving operators. In general, since these operators are confined to the branes, 
volume suppression will render the associated corrections small . 2 For a numerical estimate 
of their effects on the KK spectrum and thence the Higgs potential, see Appendix F. 

Another interesting possibility is the presence of Fayet-Iliopoulos (FI) ZVterms on the 
branes. When the hypercharges of the (tree-level) massless scalar spectrum satisfy JY Y t 7 ^ 0, 
there is a quadratically-divergent one-loop diagram leading to incalculable FI terms localized 
on the orbifold fixed points [35, 40, 48, 49]. These FI terms have two effects: first, they lead 
to explicit Higgs soft masses through the .D-term potential; second, they quasi-localize all 
hypercharge-carrying states to the 4D branes through their 5D gauge kinetic terms. The 
latter (subdominant) effect modifies the 5D wave functions, changing the effective couplings 
of these states to the Higgs through their wave function overlaps. In principle, these effects 
can be included as a free, incalculable parameter in the model. In the analysis below, we 
will assume the quadratically divergent contributions are not present and will set them to 
zero. This is consistent with the low-energy degrees of freedom for the model in which the 
Higgs states are localized to the brane at y — 0 (see Sec. IV). For the model with the Higgs 
propagating in the bulk, such a quadratic divergence does arise. However, as we will show 
in Sec. V, the light Higgs doublet lies along a .D-flat direction; consequently, the FI terms 
do not contribute to its soft mass. 

III. Folded Supersymmetry 

Having reviewed the general features of SUSY breaking by boundary conditions, this 
section will explain how the folded approach yields a model of neutral naturalness. We will 
begin with a review of the original incarnation of /-SUSY, which is a 5D model with maximal 
Scherk-Schwarz twist [4], We will provide an estimate indicating that /-SUSY at maximal 
twist does not break electroweak symmetry, motivating the exploration of models with non¬ 
trivial twists. We then discuss the general features of these models, which yield a non-zero 
Higgs vev at one loop. This section provides the requisite background and motivation for 
the two detailed models that we will consider in Secs. IV and V. 

The essential idea underlying this avatar of neutral naturalness is the pairing of opposite- 


2 As discussed in Appendix F, there is a possible exception when brane-kinetic terms are allowed to be 
negative. In this case, the existence of KK states with zero norm at tree level necessitates working to 
higher order in perturbation theory. We leave a full analysis of this scenario to future work. 
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TABLE I: Transformation properties at maximal twist. Note ^(j) is a generic matter (/-)fermion, 
'L(j) is a generic matter field (/-)scalar, A M is a generic gauge boson, E is the associated scalar 
that fills out the AT = 2 vector multiplet, and A 1.2 are the two gauginos. The brackets denote the 
choice of boundary condition for each field under Z 2 and Zl 2 respectively. 

spin partner states with different Standard Model quantum numbers via discrete symmetries. 
Since 4D J\f = 1 global supersymmetry necessarily commutes with the Standard Model gauge 
groups, this pairing looks like a hard breaking of SUSY. It may, however, be successfully UV 
completed via a 5D supersymmetric embedding with SUSY-breaking boundary conditions. 
Discrete symmetries relate bulk fields charged under different gauge interactions, ensuring 
the equality of their couplings. Dimensional reduction of this theory leaves behind the 
desired spectrum of light states, supplemented by towers of Kaluza-Klein modes. 

A. Folded Supersymmetry at Maximal Twist 

In the minimal incarnation of /-SUSY, the bulk gauge group is SU{2) w x f/(l)y x 
SU(3) C x 577(3)/. The matter content of the MSSM is doubled; these are the /-states, 
and they carry the same charges as their MSSM partners, with the exception that the 
colored states are charged under 577(3)/ instead of SU( 3) c . A Z{ symmetry is imposed 
that exchanges the MSSM multiplets <f> for their folded partners $/. This Z 2 enforces that 
the Yukawa couplings to the Higgs are the same for matter and folded matter and that 
the SU(3)f gauge coupling is equal to a s (at the scale where Z 2 is a good symmetry). Of 
particular importance here are the /-tops and /-stops which (along with the tops and stops) 
will yield the dominant contributions to the Higgs potential. 

At maximal twist, one can label the various fields by their eigenvalues under the Z 2 
transformations that reflect y about the orbifold fixed points at y — 0 and n R , as defined in 
Eq. (1) above. The transformation properties are given in Table I. The “folded” multiplets 
have an inverted hierarchy between the fermionic states and their superpartners, enforced 
by choosing opposite eigenvalues of Z 2 for matter versus /-matter. The resulting spectrum 
is sketched in Fig. 2. Note that the gauge sector is not folded; the zero-mode gauge bosons 
are massless and the lightest gauginos have mass 1/(2 5). 
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FIG. 2: Sketch of KK top sector spectrum at maximal twist. 


The following Yukawa couplings are allowed: 


WD6(y)yi 5) (HQU + HQ f U f ), 


( 5 ) 


where is the 5D top Yukawa coupling, and imposes that the Yukawa couplings are 
equal. In general, Scherk-Schwarz SUSY breaking leads to tree-level non-supersymmetric 
splittings within bulk multiplets, which would then feed down into one-loop soft masses 
for brane-localized states. At maximal twist, however, the one-loop Higgs potential enjoys 
an apparent accidental SUSY due to “bi-fold” protection [4] - the scalars from the folded 
sector /-top supermultiplets are degenerate with the fermions in the Standard Model top 
supermultiplets (and visa versa). This leads to a cancellation of the one-loop threshold 
corrections from Eq. (5) for the Higgs mass at each level of the KK tower. However, this 
bi-fold protection holds only for the contribution of top-stop loops to the Higgs doublet 
masses. Gauge invariance of the Yukawa couplings requires all states to be charged under a 
single SU{2)w x U(l)y, and so there can be no doubling in the gauge sector. Accordingly, 
clectroweak gauginos contribute to the Higgs mass at one loop, as in the 5D MSSM. Thus 
at maximal twist the Higgs potential accumulates one-loop contributions from the gauge- 
gaugino sector but not the stop-top sector. 

The squarks likewise receive soft masses at one loop from their gauge and Yukawa 
couplings. In particular, the (/-)stop soft masses are [4] 
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where K = 7£(3)/4 ~ 2.1 and an exact l{ symmetry on the values of the couplings has 
been assumed for the numerical evaluation. 

At two loops, these SUSY-breaking masses feed into the Higgs potential. It is reasonable 
to assume that the dominant effect of this two-loop soft-breaking of SUSY can be captured 
by including the one-loop soft masses as a shift in the pole of the propagator for the relevant 
states. 3 Note that we model this effect as a universal shift in the mass of the entire tower 
of (/-)top squarks, which is a reasonable approximation given the explicit values shown in 
Eqs. (6). Then the contribution to the Higgs mass squared parameter can be computed 
using standard 5D techniques, as described in Appendix D. The result is 


2 2 , 21 C(3) (g 2 + g ,2 /3) , a2 

m H = /I d- „, , ^ -1- A 


64 7T 4 R 2 


mr 


3 (log [(4 7 xR mi) 2 ] - 3) 


8 7T 2 


+ 


(7) 


where /j is the standard MSSM Higgsino mass, the second term is due to the one-loop 
contribution from gauge/gaugino loops, and the third term comes from our approximate 
evaluation of the leading two-loop effects. We denote by rrq the universal one-loop stop//- 
stop soft mass in the approximation that all the stops and /-stops receive identical radiative 
corrections to their masses. 

Using Eq. (7), it is straightforward to see that achieving m 2 H < 0 requires nz~ > 0.04/i? 2 
when fj, = 0 (which would be excluded by chargino limits if the Higgs is brane-localized, 
but is allowed if the Higgs propagates in the bulk and the Higgsino is rendered massive by 
Scherk-Schwarz boundary conditions); rn 2 needs to be strictly larger than this for /i / 0. 
Combined with Eqs. (6), this provides evidence that the model as originally constructed 
does not break electroweak symmetry. Specifically, the one-loop stop//-stop mass would 
need to be a factor of ~ 3 larger before this two-loop effect would overwhelm the one-loop 
soft mass from the gauge sector and drive electroweak symmetry breaking. 

We note that our approximate two-loop contribution in Eq. (7) differs from the two-loop 
estimate in [4], which was used to argue for the viability of electroweak symmetry breaking 
at maximal twist. The estimate in [4] was performed by running the Higgs mass from the 
(/-)t°P sector states starting at the scale 1 /R, yielding a dependence on log (Ruif). However, 
the calculation involving the full sum over KK modes depends on log(4 tt R my), which proves 
critical for overcoming the positive gauge contribution. This can be understood from the 
five-dimensional perspective - the argument of the logarithm should be a function of 2 7r R, 
rather than A, because loops must propagate a distance 2 7 t R from one boundary to the 


3 This is in close analogy with “daisy resummation” in finite temperature field theory (see e.g. [50] for a nice 
review), except here we are compactifying a spatial dimension from 5D to 4D as opposed to compactifying 
the time direction from 4D to 3D. 
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other and back again in order to see that SUSY is broken. Note that the finite terms in 
Eq. (7) also play a numerically non-trivial role. 

There are a couple of important caveats to this claim. It is possible that additional 
corrections from the full two-loop calculation could yield the proper vacuum. Furthermore, 
this calculation assumes that the incalculable, brane-localized Kahler potentials for the top 
and /-top superfields are negligible. Nevertheless, we take this as motivation to extend the 
model in order to achieve electroweak symmetry breaking at one loop. 


B. Folded Supersymmetry at Arbitrary Twist 


The choice of a = 1/2 made the top quark and /-stop zero modes degenerate at tree level, 
thereby endowing /-SUSY with bi-fold protection (see Fig. 2). By moving a away from 1/2, 
the degeneracy is broken, and with it, the bi-fold protection of the Higgs potential. 

The twist can be formulated as follows. There are two possible actions for 5D orbifolds: 
reflection Z about y = 0, and translation T by 2n R. These can be identified with the 
Z reflection about y = 0 and the Z' = T Z reflection about y = n R defined in Eq. (1). 
Discussing the twist in the language of the translation is more natural from the perspective 
of action on superfields which are not eigenstates of Z' 2 . 

For the matter fields, the Z' 2 acts as 

^ ( t* ) ( »> = ( 4 t,) ^ R ~y) = T ■ ■ (£) <»>■ < 8 > 


where 



afT2 (— l)f for SU(2) r doublet 

1 for SU (2) f { singlets 


( 9 ) 


and / = 0 (1) for matter (folded) fields. The gauginos transform as do the matter fields in 
Eq. (9), while the A M and S transform as in Table I. The transformation properties of bulk 
Higgses are given in Appendix B. 

After electroweak symmetry breaking, the combination of the non-maximal twist and the 
top Yukawa on the y — 0 brane yields kinked KK wave functions for the stops: 
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FIG. 3: Sketch of KK top sector spectrum at arbitrary twist a. The arrows indicate how the 
masses of the states change as a deviates from 1/2. 


with { y } = y mod 2itR, such that {y} is restricted to the fundamental domain. Here 
k n ± = n + (1/7r) tan _1 (± tt Ry t v), n e Z, y t is the 4D top Yukawa coupling, and 
v ~ 174 GeV is the Higgs vev. The wave functions for the /-stops are given by taking 
Eq. (10) and replacing a —> a — 1/2. These wave functions can then be used to compute the 
Coleman-Weinberg potential from top/stop and /-top//-stop loops. The resulting tree-level 
masses for the KK towers of (s)tops and /-(s)tops are 


(n) 

m r = 


k n + a 
R 


in) k n T 1/2 
mt f R 


(n) 

m'r = 

t f 


k n + 1/2 — ol 
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This spectrum is sketched in Fig. 3. 

Since it acts differently on SU (2 ) R doublets and singlets, the non-trivial twist also impacts 
the gaugino masses. For every vector supermultiplet, there is a zero-mode Majorana gaugino 
of mass = a/R, plus towers of Majorana gauginos A*™) of masses rn^ 1 ' 1 = \n + a\/R. 

These towers arrange themselves into Dirac gauginos at maximal twist. The lightest gauginos 
consist of the usual Majorana bino, wino, and gluino, plus an /-gluino, all with (tree-level) 
mass a/R. 

In contrast to maximal twist, these a-dependent, SUSY-breaking mass splittings allow for 
clectroweak symmetry breaking at one loop. However, this comes at a price. While /-SUSY 
at maximal twist had the virtue of lifting all new colored states to ~ and preserving 
naturalness through light /-stops, /-SUSY at non-maximal twist features colored states at 
the scale ~ ^ with a < 1/2. For a given radius of compactihcation, this brings colored states 
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closer to the weak scale, where strong LHC limits enter to provide relevant constraints. As 
we will see in the following sections, it will not be possible to achieve the measured Higgs 
vev and mass without falling afoul of existing phenomenological constraints, and so some 
additional ingredients are required. Before proceeding to the complete models, we will first 
make a few generic comments about /-SUSY phenomenology. 

C. Folded Phenomenology 

Folded SUSY contains all of the phenomenology of the MSSM, along with a very rich 
set of observables from the folded sector. Given the choice of requiring calculable SUSY 
breaking from extra dimensions, we will demonstrate below that the dominant constraints 
on the parameter space will be due to standard MSSM-like signatures: Higgs mixing and 
heavy Higgs production for the brane Higgs model (see Sec. IV) and gluino-squark limits for 
the bulk Higgs model (see Sec. V). Generically, twisted /-SUSY models predict an MSSM- 
like discovery immediately followed by the observation of folded-sector states. The purpose 
of this section is to briefly review some of these observables, along with some additional 
complications that are required to make /-SUSY compatible with cosmology. 

Recall that the folded states couple to the Standard Model via the electroweak gauge 
group along with their Yukawa couplings to the Higgs (this is only relevant for the 3 rd 
generation) - thus folded states may be produced via either electroweak processes or the 
Higgs portal. Since there are no light fermions in the folded sector, the folded phenomenology 
is very different from that of the visible sector. 

The phenomenology of /-squark production at hadron colliders has been studied in detail 
in [4, 5, 18, 21, 22]; here we briefly summarize the relevant physics and make several new 
observations. When pairs of /-squarks are produced at the LHC, the lack of light fermions 
implies that the color string (with associated confinement scale A f Aqcd) connecting 
these states cannot be broken by pair production from the vacuum. The result is that the 
/-squark pairs rapidly radiate /-glueballs and soft photons before eventually annihilating. 
The dominant annihilation products depend sensitively on the properties of the /-squark 
pair. Neutral-current production of light-flavor /-squarks typically ends in annihilation to 
twin glueballs, while charged-current production typically ends in annihilation to W ±r y. 

The situation is somewhat different for /-stops and /-sbottoms. If the /-stops are 
heavier than the /-sbottoms, they typically undergo /3-decay before annihilating, and the 
phenomenology is similar to that of light-flavor /-squarks. However, if the /-stops are lighter, 
then the /-sbottoms undergo /3-decay into the /-stops, and annihilation may proceed into 
final states involving Higgs bosons. Although the spectrum of /-squarks is nominally fixed 
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in terms of a and R , the effects of incalculable brane-localized wave function renormalization 
factors can yield a re-shuffling of the /-squark spectrum, so that all possible orderings should 
be considered. Thus the generic signatures of /-squark pair production are resonances in 
VV' channels where V, V' = W ± , Z°, 7 , h, gf, depending on the identity of the /-squark 
pair. Note that the timescale for this annihilation is prompt at the LHC. Limits have been 
estimated for these channels, leading to a bound of rn Cif > 300 — 500 GeV from W ± 7 
resonance searches [5]. Although no current bound has been placed using W ± h and hh 
resonance searches, these should provide interesting limits as well. Decays entirely into /- 
glueball states followed by displaced decay to light Standard Model states are also likely to 
prove promising [ 22 ], 

There is also phenomenology associated with the /-sleptons. In particular, the lightest 
/-slepton is stabilized by /-lepton number and, in some cases, .R-parity conservation. If it 
is electromagnetically charged, there are strong cosmological constraints. A simple way to 
avoid such constraints is to add additional terms to the superpotential that violate these 
symmetries. For example, including yfj Lf Qi Dj in the superpotential on the y — 0 brane 
would allow the /-sleptons to decay to di-jets. The Yukawa coupling is included so as to 
avoid spoiling the MFV assumption for the MSSM states. As the /-slepton must only decay 
before BBN, the residual flavor violation from the /-lepton sector can be made negligible. 
Despite the .R-parity violation, the lightest neutralino remains stable, as all lighter fermions 
carry either Standard Model baryon or lepton number, which are unbroken by the addition 
of this operator. 

Alternatively, another version of /-SUSY could be constructed where the l{ symmetry 
only applies to colored states, such that there is no doubling of bulk lepton supermultiplets. 
In this case, the zero-mode spectrum remains anomaly-free despite the elimination of folded 
lepton supermultiplets as all /-leptons are made massive from boundary conditions even 
when the lepton sector is folded. The only possible consequence of the missing matter is the 
introduction of bulk Wess-Zumino consistency terms [49], which have no significant impact 
on the phenomenology discussed here. 

IV. Higgs on the Brane 

Here we consider the field configuration of the original model of /-SUSY: the matter and 
gauge fields are taken to propagate in the bulk, so that their zero-mode superpartners are 
all lifted by Scherk-Schwarz SUSY breaking, while the Higgs doublet chiral superfields H Ut d 
are confined to the y = 0 brane. This setup is sketched in Fig. 4. As in the MSSM, a //-term 
is required to lift the Higgsinos above current bounds, leading to an irreducible, tree-level 
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H u ^d 

SU(2) W x SU(l) Y x SU{ 3) c x 5C/(3)/ 
{Q, U, D, L, E) . 

(( 3 /, Uf , Df, Lf , _ 


y = 0 y = n R 

FIG. 4: A sketch of the geography for the brane Higgs model. The standard MSSM two Higgs 
doublets H U) d live on the y = 0 brane, while the gauge groups, matter, and /-matter fields are in 
the bulk. 

tension between a large Higgsino mass and a weak-scale Higgs vev. This is in addition to 
the challenges of electroweak symmetry breaking which, as discussed in Sec. Ill A above, 
motivate extending the model to arbitrary twist parameter a. 

However, in the absence of additional ingredients, the regions of parameter space that 
break electroweak symmetry are excluded by constraints on the Higgs sector. Higgs coupling 
predictions run afoul of precision measurements at the LHC, and the extra MSSM Higgs 
states are excluded by direct searches. The heart of the issue is that the leading contribution 
to the term, given in Eq. (13), arises at one loop via the combination of the brane /i term 
and the 0(a/R ) effective A-terrn induced by the Scherk-Schwarz twist. Since /j must itself 
be 0(^^) for viable electroweak symmetry breaking, this implies b M is suppressed by an 
amount ~ ytj 4vr relative to other masses-squared in the Higgs sector. Similarly, the one-loop 
contribution to the down-type Higgs mass-squared is proportional to the electroweak gauge 
couplings and suppressed relative to the up-type Higgs mass-squared by an amount ~ Vt /g 2 ■ 
One minimal solution for lifting the second Higgs doublet without introducing new 
sources of SUSY breaking is to introduce the leading “BMSSM” operator which appears 
at dimension-5 in the superpotential [27, 51]. Supersymmetry is still broken strictly by 
boundary conditions; the new contribution acts as an effective b M term after electroweak 
symmetry breaking. The BMSSM operator also modifies the quartic of the light Higgs, but 
in most of the viable parameter space this is a small correction - the dominant correction 
to the tree-level prediction for the Higgs mass comes from the radiative quartic coupling. 
The complete viable model is given by the brane superpotential 

w = 6(y) (/» HuH d -^ !//„ H d f + yf H U (QU + Q,U,)j , 


( 12 ) 
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from which we calculate the Higgs potential in Eq. (13), including the leading radiative 
corrections. The full one-loop Coleman-Weinberg contribution from the (/-)top sector has 
been included, while the one-loop contributions to the Higgs soft mass from the SU(2) W x 
U (l)y gauge bosons and the leading i^-term are computed at fixed-order. These varying 
levels of approximation can be justified as follows: the presence of the Higgs vev on the 
y = 0 brane breaks both the 5D Lorentz invariance as well as the brane interchange parity. 
In the 4D theory, this leads to non-conservation of both KK number and KK parity, with the 
result that the number of sums over KK modes grows with the number of external Higgs legs. 
Therefore, it is important to include the full resummed Coleman-Weinberg potential for any 
contributions that are not already present at tree level. There is a tree-level contribution to 
the quartic from gauge couplings, and to H u — Hd mixing from the BMSSM terms; therefore, 
we only include the fixed-order contributions from these sources. The resulting potential is 

9 Re[Li 5 (e 4i ^ a+ k + Li 5 (e 4i ^ Q -^) - 2Li 5 (e 4 ^ 7 )] 
v= 512 7T 6 R 4 

+ if + rLXb 2 Re [ c(3) ~ Li 3( e2i ™)]) (kl 2 + kl 2 ) 

3 y/fi Im[Li 2 (e 4i7r ")] 

16 7T 3 R 

+ 9 + 8 9 (kl - |w| 2 ) 2 -2^:(|w | 2 + |w| 2 )(ww + h.c.), (13) 

with (H Ut d) = v Ut d and 7 = ^ tan _ 1 ( 7 r Ry t v u ). Note that in the plots, we have also included 
the O (1/M 2 ) contributions to the potential, which are omitted here for brevity. The effects 
of these terms on the numerical results are minor. 

Electroweak symmetry breaking compatible with existing constraints can now be realized. 
In Fig. 5 we present a slice of the parameter space in the M versus /i plane, fixing a = 1/3. 
The left panel shows contours of 1 /R, fine-tuning , 4 and the physical Higgs boson mass, while 
the right panel shows contours of the pseudo-scalar Higgs mass m^, and tan /3 = v u /v d - The 
region excluded by constraints from Higgs coupling measurements is shaded purple. The 
region with 2/j, > \M\ is shaded in grey, in order to demarcate the region where the BMSSM 
approximation is expected to break down. Given that the BMSSM term only serves as an 
effective parameterization of possible extensions of the Higgs sector, there should in principle 
be additional light states that contribute to the Higgs potential in this gray region. Since 
the tuning is larger here, we do not pursue the details of EWSB for these parameters. 


(v u v d + h.c.) 


4 Fine-tuning is computed using the standard Barbieri-Giudice measure A = maxi(dlogv 2 /dlog A'j), where 
x = {y,yt,g,a,M} [ 52 ]. 
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- 7!l 'i, [GeV] — Tuning l/R [TeV] - mu [GeV] — tua [GeV] . tan/3 




n [GeV] n [GeV] 

FIG. 5: A slice of the brane Higgs model parameter space in the M versus fi plane for a = 1/3. 
The left panel shows contours of fine-tuning [dashed red], 1/1? [dotted orange], and the SM-like 
Higgs boson mass [solid blue]. The right panel shows contours of tan/3 [dotted orange], and the 
pseudo-scalar Higgs mass rriA [dashed green]. The purple region is excluded by Higgs coupling 
measurements and the gray shaded region denotes 2/x > |M|, where the BMSSM approximation 
begins to break down. The minimum allowed tuning compatible with a 125 GeV Higgs mass is 
~ 0 . 2 %. 

The dominant limits in this model are due to Higgs mixing constraints, with all parameter 
space in the purple region excluded by Higgs coupling measurements. These constraints 
exclude even larger values of than they would in the MSSM. Although the BMSSM 
operator in Eq. (12) allows the pseudoscalar mass rriA to be raised, it does so in a way that 
slows the approach to the decoupling limit relative to the MSSM. Increasingly stringent 
limits on Higgs coupling deviations, largely driven by ATLAS Higgs coupling measurements, 
translate into a substantial bound on the scale of additional Higgs states in this model. To 
characterize these bounds, we perform a combined fit to ATLAS, CMS, D0, and CDF Higgs 
measurements using the package HiggsSignals v. 1.3.2 [53], parameterizing deviations in 
terms of the tree-level mixing angles in a Type 2 2HDM. Here we include the effects on 
loop-level couplings from tree-level mixing but no additional loop contributions from new 
states. We translate these bounds into a 95% exclusion in the M versus fi plane, where the 
tree-level mixing angles are computed numerically at each point in parameter space. 

In principle, the parameter space is also constrained by direct search limits on heavy scalar 
and pseudoscalar Higgs states, which extend to ~ 1 TeV at large tan/3 [54, 55]. However, 
the viable parameter space for EWSB falls at moderate values of tan/3 with uia >1.5 TeV 
already. Constraints on MSSM superpartners could also be relelvant, but given the large 
l/R > 25 TeV, these constraints are subdominant. The result is that /-SUSY with brane- 
localized Higgses is tuned at around the 0.2% level. In the least-tuned regions of parameter 
space compatible with Higgs coupling measurements, the /-stop masses are on the order of 
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Uu Ud,e 

SU{2) W x SU(1) Y x SU(3) C x SU{3) f 
(Q, U, D, L, E). 

(<2/> U E D fi L fi E f). 

H a ,b 


y = 0 y = it R 

FIG. 6: A sketch of the geography for the bulk Higgs model. The Higgs doublets Ha,b, gauge 
groups, matter, and /-matter fields are all bulk fields. The up-type (down-type) Yukawa couplings 
are on the y = 0 (y = ir R) brane. 

8 TeV. As such, the first signatures would arise in the form of Higgs coupling deviations. 
The heavy MSSM-like states and the more exotic folded-sector phenomenology will be of 
order 8 TeV and 4 TeV, and will be out of reach of 13 TeV LHC searches. 

While the parameter space of the brane Higgs model is strongly constrained by Higgs 
couplings, /-SUSY with a twist is sufficient to put superpartners out of reach of current LHC 
searches. Of course, additional sources of brane-localized SUSY breaking would allow the 
new states in the Higgs sector to be decoupled without substantially increasing electroweak 
fine-tuning, albeit at the cost of reduced predictivity. 

V. Higgs in the Bulk 

The experimental tensions discussed in the previous section motivate exploring models 
where the Higgs propagates in the bulk, since the choice of boundary conditions that yields 
a light scalar Higgs is compatible with one that lifts the Higgsinos and other Higgs scalars to 
scales ~ 1 /R. As we will discuss below, a pair of bulk hypermultiplets Ha,b will be required. 
A sketch of the geography for this model is given in Fig. 6. Note that in this setup, the 
up-type Yukawa is confined to the y — 0 brane, while the down-type Yukawa couplings are 
located on the n R brane (and involve H^, e.g. W D H^Q D). 

Following [34, 56], a natural choice would be to attempt to fold the 5D MSSM with a 
single bulk Higgs hypermultiplet. In fact, the MSSM version of that model predicts a Higgs 
boson with mass consistent with the now-observed value [34], However, the corresponding 
prediction for the radius of compactihcation is 1/R ~ 700 GeV, with a gluino around 
350 GeV. One could hope that folding this model at arbitrary twist would yield a partial 
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reduction of the negative, ^-dependent Higgs soft mass, thereby raising 1/R. However, this 
fails because at non-trivial a, the zero-mode Higgs doublet receives a positive KK mass at 
tree level. 

One viable extension is to introduce a second bulk Higgs doublet. An additional SU(2) B 
symmetry rotates the two Higgs hypermultiplets Ha,b into each other. The twisted boundary 
conditions can be modified to take advantage of this extra freedom in order to yield a single 
massless Higgs doublet at low energies. Specifically, the action of the Z 2 at y — 0 is taken 
to be 
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while the 2 n R translation T is given by 
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Note that we have implicitly set the SU(2) B twist angle equal to the SU(2) B twist angle. 
Whether these angles are free parameters or instead determined by dynamics depends on 
physics above the cutoff. If they were unequal, then the light Higgs boson would receive an 
0(1/R 2 ) tree-level positive mass squared proportional to the misalignment between these 
parameters. 

This choice of boundary conditions leads to a single light Higgs scalar h ^ and heavier 
Higgs scalars and higgsinos with masses ~ a/R] the KK decomposition is given in Eq. (B6). 
One of the consequences of this configuration is the vanishing of the tree-level D-terrri quartic 
for A°d Therefore, in contrast to the brane Higgs model, an additional source of quartic is 
required in order to lift the Higgs boson mass to the observed value. We choose to do this 
by adding a brane-localized singlet S and an NMSSM-like coupling between the singlet and 
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the Higgs bosons: 

W = 8(y) (\®SH a H% + ^ S 2 + H A (QU + Q f U f )^j . (16) 

Note that since the singlet has a large supersymmetric mass, it will not get an appreciable 
vev; furthermore, singlet-Higgs mixing will remain sufficiently small so as not to affect Higgs 
coupling measurements. 

For the parameter space of interest, the singlet-Higgs coupling will be 0(1) to accommo¬ 
date rrih = 125 GeV. In particular, the Higgs mass-squared is entirely determined by one- 
loop contributions, predominantly from the singlet and top Yukawa couplings (neglecting 
the small tree-level contribution from the singlet vev). Therefore, it is important to include 
the leading radiative contributions to the Higgs potential from these couplings. Following a 
similar justification to that given in Sec. IV above, the full Coleman-Weinberg potential from 
the (/-)top states is included. The one-loop contribution to the Higgs mass parameter from 
loops involving the singlet is evaluated at fixed order, as there is a tree-level contribution to 
the quartic proportional to A^. Note that there is no longer a tree-level quartic from the 
SU (2) vc x U (l)y gauge bosons. However, the Higgs vev no longer breaks KK number/parity 
(although the top Yukawa coupling does), removing the proliferation of sums over KK modes 
in gauge/gaugino loops. Thus we need only include the fixed-order contribution to the Higgs 
mass. In analogy with Eq. (13), the potential for the zero-mode Higgs is 

V = (A v 2 — M s vs ) 2 + 2 A 2 Vg v 2 

9 Re[Li 5 (e 4i ^ Q+7 )) + U 5 (e 4i ^ a ~0) - 2Li 5 (e 4i7r7 )] 

512 vr 6 /? 4 

+ (ifwrfr Re [C(3) - Li 3 (e 2l7r ")] + m^ Singlet ) u 2 , (17) 

where = v, ( S ) = vs, 7 = \ tan -1 (tt Ry t v), m 2 H Singlet is the soft mass contribution 

from the singlet sector (see Eq. (E5) in the appendix for the integral form), and A = 
A( 5 )/(4 ttR) is the rescaled singlet-Higgs coupling. Note that S acquires a tadpole that leads 
to a small vev vs = (Mgv 2 A)/(Mj + 2r? 2 A 2 ), which is included in the numerical results. 

Now that we have the tools to explore the vacuum structure of the model, we may 
work out the dominant constraint on the parameter space. Recall that in this model the 
gauginos, squarks, and Higgsinos are all bulk fields. Given the choice of boundary conditions, 
the zero modes for all of these states receive a/R- sized masses from boundary conditions. 
Naively, this implies that the gluinos, squarks, and neutralinos will all be degenerate, while 
constraints on gluino-squark-neutralino simplified models will yield the strongest bounds. To 
determine precise bounds, we must consider both the dominant gluino decay channels and 
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the level of degeneracy between the various states. The 0(a/R ) estimate for the neutralino 
mass is modified by mixing due to the vev of the Higgs (along with the presence of the 
S'-fcrmion, whose mass is similar to a/R in the viable parameter space). The neutralino 
and chargino wave functions can be computed explicitly including the discontinuities that 
result from the Higgs and singlet vevs. We find that the neutralino mass is degenerate with 
the colored superpartners with a splitting of at most 0(20%). 5 

In addition to calculable contributions to the masses, there are incalculable contributions 
from brane-localized Kahler terms that can shift the mass spectrum as discussed in detail 
in Appendix F. This can change the qualitative features of the LHC signatures in various 
ways. While we expect the splitting between the gluino and lightest neutralino to be robust, 
for extreme cases this degeneracy could be spoiled. Another interesting variation occurs by 
changing the identity of the lightest squark. For example, if the stop becomes lighter than 
the light-flavor squarks, the gluino decay proceeds through one or two off-shell tops. Gluino- 
squark associated production then involves a decay to a nearly degenerate neutralino, two 
tops, and light-flavor jets. This specific signature has not been searched for directly at the 
LHC and presents an interesting opportunity for a new analysis. 

While this somewhat-exotic signature is possible, for a wide range of the parameter space 
the dominant channel will be gluino and squark production with decays into light-flavor 
quarks and nearly-degenerate neutralinos. Hence, jets + MET searches at the LHC are 
relevant [57, 58]. An exclusion plot in the gluino/squark versus neutralino plane in the 
associated gluino-squark production channel for this model has been given in [57]. While 
this result does not include all the possible production modes in the simplified model, it does 
provide constraints on the largest production channel and therefore provides a reasonable 
approximation of the full limit. In particular, this result captures the degree of degradation 
to the limits that results from degeneracy of the spectrum. Therefore, we use this published 
limit as an estimate of the exclusion in our results below. 

Two slices of the viable parameter space are given in Fig. 7: the left [right] panel shows A 
versus a [Ms R]. The red region is excluded from jets + MET searches at the LHC [57, 58], 
and in the grey region there is no electroweak symmetry breaking. Contours of the SM- 
like Higgs boson mass are shown in blue; the black contours denote 1/-R; and contours of 
fine-tuning 6 are given in brown. We can see that while the LHC has probed this parameter 
space, there is still a wide region where the tuning is mild (A < 10). Jets + MET searches 

5 Note that the connection between degenerate spectra and 5D models was made previously in [33], although 
there the Higgs was brane-localized so that the level of compression, which depended on the value of the 
/x parameter, was less robust. 

6 The fine-tuning measure is described as in the previous section with X = {y t , g, A, a, Ms}. 
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FIG. 7: Brane Higgs parameter space in the singlet-Higgs coupling A versus the twist a plane 
[left] and the A versus the singlet mass parameter Ms plane [right]. The blue contours denote 
the mass of the Standard Model-like Higgs, while the black contours denote the inverse radius of 
compactification l/R. Contours of fine-tuning are shown in brown. The red shading indicates the 
exclusion coming from gluino mass limits relevant to the compressed parameter space of the model, 
while the grey shading indicates the region where EWSB does not occur. 

at LHC13 will provide a strong test of the minimally timed region of parameter space, 
and the observation of /-squarks would follow a discovery in this channel. While the level 
of fine-tuning here is similar to that found in other effective theories of “Natural SUSY” 
with low cutoffs, it is worth emphasizing that /-SUSY provides a fully calculable model of 
SUSY breaking. In comparison to other such models, e.g. low-scale gauge mediation, the 
fine-tuning in /-SUSY is quite mild. 


VI. Conclusions 

In this paper we have constructed the first calculable models of folded supersymmetry 
with successful electroweak symmetry breaking at one loop. To do so, we have generalized 
the original /-SUSY framework to incorporate non-trivial Scherk-Schwarz twist. We have 
also exploited some of the geographical freedom allowed by 5D theories, considering models 
in which the Higgs doublets are brane- or bulk-localized. 

In brane-localized models, radiative corrections from bulk states naturally accommodate 
a 125 GeV Higgs. However, the smallness of certain soft-mass parameters in the Higgs 
sector leads to additional scalar and pseudoscalar Higgs bosons excluded by precision Higgs 
measurements. This may be addressed by deformations of the theory (such as the dimension- 
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5 BMSSM operator) that raise the mass of the pseudoscalar Higgs, but even so Higgs coupling 
measurements provide the strongest constraint on the parameter space and push tuning of 
the weak scale to the sub-percent level. 

In models where the Higgs doublets propagate in the 5D bulk, constraints from the 
Higgs sector are mitigated because the additional states are lifted by boundary conditions. 
However, the viable model introduced above has no D-term potential, requiring the 
introduction of new tree-level quartic couplings a la the NMSSM. Although light folded 
states preserve the naturalness of the weak scale, the region of parameter space leading 
to electroweak symmetry breaking brings colored states into the reach of direct searches. 
The strongest bounds on these bulk Higgs models come from the gluino-squark-neutralino 
simplified model. In light of current limits in the jets + MET channel, these models are 
tuned at better than ~ 10%. 

As a whole, this work highlights both the perils and the opportunities inherent in attempts 
at building concrete models of neutral naturalness. Accommodating electroweak symmetry 
breaking and the observed Higgs mass pushes calculable realizations of /-SUSY away from 
the limit in which colored partner states are fully decoupled. Rather, it favors regions 
with an admixture of conventional and folded phenomenology. The strongest bounds may 
come from more conventional channels, including precision Higgs measurements and colored 
sparticlc searches. 

Our work also highlights the vast range of possibilities contained within the /-SUSY 
paradigm. There are a variety of worthwhile additional directions in which /-SUSY may be 
deformed while preserving the essential mechanism leading to neutral naturalness: 

• Further exploiting the geographical freedom in 5D models. The original /-SUSY model 
featured brane-localized Higgses and bulk-localized matter fields. In this work we 
have considered the implications of bulk-localized fields, but in general there are a 
variety of possible configurations of bulk- and brane-localized fields that may have 
novel consequences. 

• Exploring the opportunities provided by anomaly cancellation. Crucially, the boundary 
conditions of /-SUSY leave no light folded fermions, so that zero-mode anomaly 
cancellation is insensitive to the precise choice of bulk folded states. Consistency 
of the 5D theory still may require the introduction of bulk Chern-Simons terms, but 
these do not significantly alter the folded mechanism. We have pointed out that this 
may be used to evade bounds on folded sleptons, but may be exploited more generally 
to produce more minimal models along these lines. 

• Considering additional sources of SUSY breaking. In this work we have restricted 
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ourselves to SUSY breaking by boundary conditions in order to preserve the 
calculability of the Scherk-Schwarz twist. However, this restriction leads directly to 
many of the dominant constraints in our models, which may be ameliorated by the 
introduction of additional sources of SUSY breaking. In particular, this may imbue 
the original folded proposal with viable clcctroweak symmetry breaking and weaken 
bounds from conventional searches for SUSY partner states. 

In the absence of new physics signals, the question of the smallness of the weak scale 
becomes increasingly pressing. The exploration of novel theories that explain the lightness 
of the Higgs mass continues, and neutral naturalness represents a promising paradigm that 
accommodates null results but nonetheless yields experimentally observable consequences. 
As we learn more from the now-ongoing run of the 13 TeV LHC, data will only sharpen the 
question and encourage further exploration. This work has demonstrated that if we live in 
a folded world, then it is likely that we will know about it very soon. 
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A. Boundary Conditions 

The starting point for all of our calculations will be to solve for the boundary conditions 
of the bulk fields. At arbitrary twist angle a, doublets of SU(2) R obey the following 
transformation properties: 

(**) ) M < A1 > 

where / = 0 (1) for the MSSM (folded) fields, and a* are the Pauli matrices. Any generic 
SU(2) r twist can be brought to the form in Eq. (Al) through an SU(2)r transformation. 

For fields that do not transform under SU(2)r, we define an “effective” Scherk-Schwarz 
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twist. For example, the bulk matter fermions have the following transformation properties: 

(j)(» + 2 ' JI ) = t- 1 ) / (j,)(»); (£)(-») = "* (£)<»)• (A2 > 

where (ij), ij) c ) are the fermionic components of a bulk hypermultiplet. Since e~ l7r ^ a2 = (—lp , 
the MSSM (folded) fermions have an effective Scherk-Schwarz twist a = 0 (1/2). 

In order to simplify calculations involving both MSSM and folded bulk holds, the result 
can be obtained by assuming there is only one SU (2 )r doublet with twist a. Then the 
contributions from the additional states is captured by replacing a with the appropriate 
effective Scherk-Schwarz twist listed in Table II. 



MSSM 

Folded 

Fields 


(A?,*A£),(Q,Q c t) 

( Qf,Q }) 

«/,<??) 

Effective SS-twist 

0 

a 

1 

2 

a ~\ 


TABLE II: The effective Scherk-Schwarz twist for different 5D bulk states. The results for a given 
field can be extrapolated from those for SU{2)r doublets by replacing a with the appropriate 
effective Scherk-Schwarz twist. Note that for fields with zero effective SS-twist, an extra 1/V2 
factor must be included to appropriately normalize the 4D kinetic term. 
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B. Wave functions 

Here we perform the KK expansions for all bulk fields. At non-maximal twist, squarks 
and gauginos mix with their SU(2) R partners. Their KK expansions take the form of a sum 
over SU(2)r doublet wave functions, ^Sf n (y): 


E Q {n \^)^n(y) 


E X ‘n)( X “) «’»(»)• 


Here ^ n (y) needs to satisfy the boundary conditions in Eq. (Al), and is given by 


*n(y) = 


(»+<*) y 


'2-rR I sin ( Mp 


where n is an integer. The mass for each mode is m n = \n + a\/R. The fermionic wave 
functions and their folded counterparts can be obtained via Table II. 

In general, Scherk-Schwarz symmetry breaking boundary conditions can involve symme¬ 
tries other than SU (2 ) R , this can lead to more complicated wave functions. For example, for 
our model where the Higgs fields are in the bulk, the Scherk-Schwarz twist can be embedded 
within any element of the SU{2) R x SU(2) H symmetry. One can fix an SU{2) R and SU{2) H 
basis such that the transformations for the Higgsinos are 


n A n A 

H b H c b 


(9 + 2 = H - A (9); 

\ Hb R b 


' ha hj 

H b HI 


' Ha H c a 

H b H% 


^3 (y)- 


and for the Higgses 


H a Hr 


ff c t u c t 

n A n B 


A B 1 _|_ 2 7 T R) — e _27riaCT 2 I A B ] e 2ir i {3 (?2 


H a Hr 


ff c t ff c t 
n A n B 


(h a h b \ 

{Hi Hi) 

The KK expansion for the Higgsinos are 


IIA H, 


^ (s) ' 


E H^ n \xA*n(y)\ 


E H c(n) (z'O'Ms/), 


where the wave functions ^ n {y) are the same as those for the squarks given above in Eq. (B2) 



with a —> /3. The ( H^ n \ H c ( n ^ pair up to form a tower of Dirac fermions with masses |n + 
f3\/R. For the Higgs bosons, the KK expansion mixes all four scalars. Ignoring clectroweak 
symmetry breaking, we find 


H a H b 

jj c t u c t 

n A n B 


£ 


/i<”V)M») +if ( ”V) «■,(») 


where the matrices h n (y) and ti n (y) are given by 


( cos ( 


h n (y) = 


R 


V R 

•sin( (w+a ~^ )y 

((n+a+0) y 


sin ( b±Zz®v \ \ 


cos ( (n+a ~^ )y 


/ 


cos 


H„(t/) = 


V R 


VMi l _ sln ( fr+a+ftv j _ CQS ^ ^ 


sin \ 

((n+a+0)v\ 

/ 


(B6) 


(B7) 


Before clectroweak symmetry breaking, the masses for the h^ and modes are n + 
a — /3\/R and \n + a + f3\/R, respectively. Fixing a — (3 results in a massless h ^ which is 
desirable as a starting point for models of natural electroweak symmetry breaking. From 
the 4D perspective, h ^ is the SM-like Higgs held that develops an clectroweak symmetry 
breaking vev v. This leads to additional w-dependent contributions to the mass spectrum 
not explicitly shown here. Note that because Ha and H'^ contribute equally to h^°\ its 
potential will not receive contributions from the gauge D-terms. 


C. Boundary Fields and Discontinuities 

When the 5D Lagrangian contains interaction terms localized on the boundaries, the 5D 
wave functions can develop kinks and discontinuities. One such important interaction comes 
from the Higgs; a non-zero vev for a 4D Higgs held (either a boundary or KK state) ( H) = v 
induces a brane-localized mass. Consider the superpotential for the top sector at finite v, 

W D -6(y)(2nRy t v) U T + U c d 5 U + T c d 5 T, (Cl) 

where y t is the 4D top-yukawa coupling. The same brane-localized interaction exists for the 
brane Higgs and the bulk Higgs scenarios (v is replaced by v u in the brane Higgs case). To 
avoid any ambiguities when integrating across the 5 th dimension due to the 5(y) function, 
we will work with the interval (—e, 2tt R — e), where e > 0 is small. A Higgs vev induces 
mixing between all four helds (t,u c ,u\t c ^), and one expects their KK expansions to be a 
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sum of vectors involving all four fields, 


ft\ 

u c 

v) 

IW 



t;(») 


(C2) 


Here T ^(y) are the 5D wave functions that satisfy the equation of motion in the domain 
y G (—e, 2 tt R — e): 


(d 5 + A^)(d 5 -A)T ± 



A 


■(2 7 rRy t v) 6 (y) 


( 0 
1 
0 

V° 


0 0 0 \ 
0 0 0 
0 0 0 
0 1 0 ) 


(C3) 


The 5(y) terms in A will force T to develop a discontinuity at y — 0. Due to the 
Scherk-Schwarz boundary conditions, the discontinuity is extended to all integer multiples 
of 2 7T R. Away from the discontinuities, the equation of motion are still solved by the 
cosine and sine functions in Eq. (B2), up to n-dependent phase shifts. In order to introduce 
a periodic discontinuity to the wave function, it is convenient to introduce the function 
{y} = y mod (2 7ri?); {y} maps the real line onto the fundamental domain [0,2 7 t R) and 
contains a finite jump from 2 7r R to 0 at every integer multiple of 2 7 r R. 

Doublets of SU(2)r that satisfy the boundary conditions and have a discontinuity at 
y — 0 can be written as 


*k(y) 


1 

V2nR 



k{y}+ay 

B 


7T k 


k{y}+ay _ ^ 


{y} = y mod 2 7r i?, (C4) 


where k is not restricted to integer values. The extra — tt k phase ensures that cr 3 ^> k (—y) = 
\E ’k(y)- Similar wave functions have been considered before in [59]. 

Note that >&*.(— y) matches the continuous version in Eq. (B2) when k is taken to be an 
integer. A sketch of the two components of the wave function \I *k(y) with a = 1/2 is shown 
in Fig. 8. The cosine component develops a kink but remains continuous since it is even 
under reflection about y — 0. On the other hand, the sine component develops a finite, 
discontinuous jump across y = 0 and is odd under reflection about y — 0. 

Taking the derivative for the sine component, the discontinuity at y — 0 gives a delta 
function. In order to satisfy the equation of motion, the 5(y) terms from 85 and A must 
cancel, which leads to an equation for k that has discrete solutions k n . The KK modes are 
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Even Odd 




FIG. 8: Sample discontinuous wave functions shown in arbitrary units as a function of y/R, 

for a = 1/2. The left (right) panel shows the shape for even (odd) boundary conditions at the 
y = 0 brane. The black, red, and blue lines denote \1/ k(y ) for \k\ = 0,1/8, and 1/4, respectively. 
The solid (dashed) lines correspond to positive (negative) k. 


found to be 


Td = 


yflVR 



k n {y}+ay _ 7Tk \ \ 
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k n {y}+ay 
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k n {y}+ay 
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— Ti k, 
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k n = n ± — tan 
7r 


-l 


( 7 vRy t v), (C5) 


where n is any integer. The operator 85 — A acting on T simply picks out the 
smooth part of the derivative; thus the masses of the KK modes are given by 
= \n + a ± tan _1 (7ri??/ f v)/tt\/R. 


D. One-loop Coleman-Weinberg Potential 

In order to evaluate the vacuum structure of the models discussed in the main text, the 
Coleman-Weinberg contribution to the Higgs potential from the top//-top sectors must be 
evaluated. At one loop [60] 

Vcw(v) = i ( -1 )^ Trlo § (<9 2 + m 2 i (v)\ 

*SDOF 


(Dl) 
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where F — 0 (1) for scalars (fermions), and v is the Higgs vev. In principle, the summation 
is over all degrees of freedom that couple to the Higgs, but here we will focus on the 
contributions proportional to y t . If the states running in the loop descend from bulk fields, a 
summation over KK modes is required. For the (/-)top sector, the v dependent KK masses 
in Eq. (C5) are analytic, and the one-loop Coleman Weinberg potential can be computed 
analytically: 


V, 


yt 


6 

If 



d 4 f 


£ 


f + a+ T) | + 77 

(f Fa + yf + F - (f+ //) 2 + £ 2 

§ + a - y | - y 

(|+a-j () 2 + F + (|-i tf+P 


(D2) 


where we have combined the contributions from the top and /-top sector together under 
the same summation over n for convenience. The potential can then be evaluated as sums 
of polylogarithms, 


9 Re [Li 5 (e 4iir ( a+ F ) + Li 5 (e 4 * ^ ^ a_7 ^) - 2 Li 5 (e 4i7r7 )' 

Vyt = 512 vr 6 IF ’ 

7 = —tan _1 (7r Ry t v), (D3) 

71 

where constant terms have been dropped. Note that at maximal twist, V yt becomes v- 
independent (at one loop) due to the bifold protection. 


E. One-loop Fixed-order Calculations 

In addition to the Coleman-Weinberg potential from the top sector, there are other 
important one-loop contributions to the Higgs potential which are discussed in the following 
subsections. 


1. Gauge soft mass 

The bulk gauginos obtain tree-level, SUSY-breaking masses from boundary conditions. 
These in turn yield a contribution to the Higgs soft mass at one loop. For convenience, we 
explicitly compute only the gaugino loops (see Fig. 9); the contribution from the bosonic 
gauge sector can be extracted by the replacement a —* 0 in the gaugino result. 

When the Higgs is brane-localized, it only couples to the J\f = 1 vector multiplet that 
is even under reflection about the y = 0 brane, such that only the diagram involving 
(A“, h ) is present. For the bulk Higgs case, there are multiple diagrams involving different 
combinations of gauginos and Higgsinos. 
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Ha,H% 



h<® 


FIG. 9: Gaugino contributions to the zero-mode (bulk) Higgs soft mass. Other contributions from 
gauge interactions and D-term quartics can be obtained from this diagram from the limit a —> 0. 
Analogous diagrams contribute to the soft mass in the case of the brane-localized Higgs. 


A straightforward diagrammatic computation reveals that the gauge soft mass in the 
brane and bulk Higgs models are the same (ignoring possible, small corrections rs_/ (/nR) 2 ): 


777 ^ 

i/, gauge 


3 g 2 + g' 2 
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3 g 2 + g' 2 
16 7r 4 R 2 
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Unlike the bulk matter contributions, the gauge soft mass is largest at maximal twist since 
there is no bi-fold protection in the gauge sector. 


2 . 


cU ct 


H u \ 


; n d 


u 


<9 5 i/H 



t 


H d 


FIG. 10: One-loop contributions to b M from the stop sector. Analogous diagrams involving the 
/-stops are also included. The gauge contribution is numerically subdominant and not shown here. 


The brane Higgs model presented in Sec. IV of the main text relies on a fi term in the 
brane-localized superpotential in order to give the Higgsinos non-zero mass. Since the tree- 
level spectrum exhibits SUSY breaking, a b M term is generated at one loop. There are two 
contributions to 6 M , one proportional to the top Yukawa and another from gauge interactions. 
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The (/-)squark contribution relies on the following interactions from the F-terrn potential: 

C D J d b x5(y) H u (^t d§u c ^ — ud^t c ^j — yyf'* Hd^ u 1 + h.c. . (E2) 

Figure 10 shows the one-loop contribution to the b^ term from the top/stop sector. Resolving 
the S(y) in these interactions requires a double summation over the KK tower - the resulting 
contribution dominates over that from gauge loops, since the latter only contains a single 
sum. Explicitly, 

_ 3 y\ n f d 4 £ T ni + n 2 + 2a f l\ 

(2W 4 ^ L [(" 1 + a ) 2 + n [("2 + a) 2 + £ 2 ] + r^ a ~2j_ 

= (E3) 

From this result, it is straightforward to see that b M vanishes at a = 0,1/2, since in both 
limits the different components of an SU(2) R doublet do not mix, i.e., a U(1)r symmetry 
is preserved. At a = 1/4, b^ is also zero as a result of an accidental cancellation between 
the top and /-top sector; the dominant contribution to 6 /t comes from the gauge sector for 
this choice of a. As discussed in Sec. IV above, the parameter space near a = 1/4 is not 
phenomenologically allowed due to large Higgs mixing. We are therefore justified in ignoring 
the gauge contribution to b M in the viable regions of parameter space. 

Naively, the convergence of this loop expansion is not particularly well-behaved as 
one includes additional insertions of the Higgs vevs. Operationally, the brane-localized 
top Yukawa comes with a S(y), yielding a proliferation of KK sums. This contribution 
should, in principle, be resummed via a Coleman-Weinberg calculation in the 2D space of 
vevs (v u , Vd). However, the parameter space that results in natural clectroweak symmetry 
breaking requires that y <C 1 /R; higher order effects involving yv u are expected to be 
subdominant. Furthermore, in the viable model presented in the main text, there is a tree- 
level contribution to an effective b M which dominates the physics of clectroweak symmetry 
breaking. Therefore, we conclude that this level of approximation is sufficient for the 
exploration of the parameter space discussed in Sec. IV. 

3. Singlet Soft Mass 

In the bulk Higgs model, a brane-localized singlet S with an NMSSM-like coupling to 
the Higgs states is included in order to reproduce the measured value of the physical Higgs 
mass. In the parameter space of interest, this coupling is large enough that the leading loop 
correction involving this interaction should be included. 
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FIG. 11: One-loop contributions to the Higgs soft mass from the singlet sector in the bulk Higgs 
scenario. 


The superpotential for this model includes 

W D S(y) SH a H%+ M ^ S 2 ^j + H c a d 5 H A + H C B d,H B , (E4) 

which can be used to construct the one-loop singlet-Higgs diagrams shown in Fig. 11. 
Explicitly, this yields 

2 A 2 f°° d££ 2 A£ 2 sin 2 (na) coth(7r£) — 2 £M S R sin(2 7r a) 

Singlet = y o M 2 R 2 + £ 2 [ COsh(2 7T £) — COs(2 7T Ct) + 

2 Mg R 2 sin 2 (2 n a) coth(7r £) + i Ms R sin(4 n a) 
cosh(2 7T £) — cos(4 n a) 

(E5) 

where we have defined A = A^ 5 -*/(4 tt R) as the effective 4D coupling. The integral is evaluated 
numerically for the results provided in the text. Since the lightest states in the singlet/Higgs 
sector are bosons (in contrast to the top/stop sector) nr^ Singlet is positive. For the parameter 
space of interest, the positive soft mass from the singlet coupling is smaller in magnitude 
than the soft mass from the top sector. 


F. Brane Localized Kinetic Terms 


The 5D theories of interest become strongly coupled at scales A not too far above 1/R. 
One way that this non-perturbative physics manifests is in the form of incalculable brane- 
localized Kahler potentials, 


/C D 5{y) Zo,Oi + 5 (y - tt R) Y. z 'oR” < F1 > 

i 3 

where O t ( O'j ) are Kahler operators that respect the local J\f = 1 supersymmetry at 
y — 0 (nR). Ignoring gauge interactions, the lowest-dimension operators are of the form 
5{y) Zij (E,- + 5(y — tt R) Z b $'■, which could, in principle, lead to sizable corrections 
to the Higgs potential. Note that, in general, these arbitrary brane-localized kinetic terms 
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lead to large rates for flavor violating observables. We will assume that the flavor structure 
of these terms is controlled by MFV' [61], such that we can neglect these constraints. 

It is also possible that large Z-factors which descend from strong dynamics could spoil 
the feature of calculability espoused above. Using naive dimensional analysis (NDA), the 
sizes of these Z-factors are estimated to be ~ 3 7 t/(2A) ~ R [62], In this section, we will 
study their impact with a focus on corrections to the Higgs potential. 

Of particular concern are the Z-factors at y — irR. The symmetry that exchanges 
matter superfields with their folded counterparts is not respected by boundary conditions 
on this brane; as a result, no relationship between matter and folded Z-factors should be 
assumed. The brane-localized Z-factors kink the wave functions for the relevant states, 
modifying the Higgs/(/-)squark wave function overlap at y — 0 and splitting the effective 
top Yukawa couplings between the MSSM and folded sectors. This could spoil the bi-fold 
protection of the KK spectrum, potentially leading to large corrections to the Higgs potential. 
Nevertheless, the underlying supersymmetry of the UV theory prevents the re-introduction 
of quadratic divergences. 

In practice, the discontinuities derived from the Z-factors lead to transcendental 
expressions for the KK spectrum, spoiling summability of the KK contributions in terms 
of simple analytic functions. To keep numerical calculations tractable, we only consider 
in detail the effects of adding a Z-factor for the U type fields at the y = tt R brane, 
corresponding to adding the following terms to the Kahler potential: 

/C D Z'u 5(y — tt R) | U' | 2 + Z' Uf 5(y - vr R) \ U' f | 2 . (F2) 


Here the prime on U' denotes the superfield whose scalar component is even around the 
y = tt R brane. In the following, we will only focus on the scalars as the fermionic wave 
functions can be extrapolated from the scalar ones by taking a —>• 0. The equations of 
motion for v!, the scalar component of U', are 


(l + Z'u 5{y — 7r R))d 2 u — d\u = 0 d 2 u c ' — <9 5 


<9.5 W 


Cf 


= 0. (F3) 


[l + Z^5(y-nR) 

The 5(y — tt R) function will introduce discontinuities in the 5D wave functions. The ansatz 
from Eq. (C4) applies, and the resulting KK expansion is given by 

, with + 


u 

u 01 1 


2 R 


= tan(7r k n ). (F4) 


2-itR 


The mass of each KK mode is given by m n = \k n + a\/R. Figure 12 shows the change in 


7 The Q c superfield is required to have the conjugate flavor charges of Q. 










36 



Z'u 


FIG. 12: The mass shift 5mo of the n = 0 KK mode in units of 1/R as a function of Z\j for a = 
0,1/8,1/4,1/2. Here we have taken AR = 5. At a = 0, dmo is identically zero by supersymmetry. 

m o as a function of Z' v . At a = 0, the mass shift is exactly zero due to the restoration of 
accidental SUSY. 

A positive Z-factor of order 3n/2A from NDA yields a 10-20% shift in the mass of the 
lightest KK mode. The mass-shifts for higher KK modes are smaller. As a result of these 
shifts, the spectrum of bulk matter fields can only be reliably computed at the 10-20% 
level. This implies an inherent uncertainty in the level of compression between the colored 
superpartners and the LSP (and in the identity of the LSP itself, in some cases). 

The Higgs couplings to the top//-top sectors also depend on the Z-factors. In fact, 
the 4D Higgs-squark and Higgs-quark coupling are not equal due to the different 5D wave 
functions at y = 0. This mismatch of couplings would naively lead to hard SUSY breaking 
and reintroduce quadratic divergences, as the loop contributions do not cancel level by level. 
However, this viewpoint mistakenly assumes a common 4D cutoff for scalars and fermions, 
which violates 5D Lorentz invariance. A better approach to study the UV behavior is to 
perform the calculations in terms of winding modes in the 5 th dimension [63], where all 
divergences come from integrals over 5D momenta. It can be shown explicitly that no 
quadratic sensitivity to the cutoff appears in the calculation of the Higgs soft mass due as 
expected to the presence of 5D SUSY. 

Next we will use the KK decomposition in Eq. (F4) to compute the modifications to the 
Higgs potential. Transcendentality of the KK spectrum implies that an analytic Coleman- 
Weinberg calculation is not possible. Hence, a fixed-order calculation will be presented to 
provide an estimate for the effects of Y-factors. 

First, consider the correction to the one-loop Higgs soft mass. The squark contributions 
come from the two sets of diagrams shown in Fig. 13. The quark contributions can be 
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FIG. 13: One-loop contributions to the Higgs soft mass from squark loops in the presence of 
brane-localized kinetic terms. The fermionic diagrams are not shown; their contributions can be 
extrapolated from the scalar ones by substituting a -» 0. 


extrapolated from these results by substituting a —> 0. The full one-loop Higgs soft mass is 
found to be 
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(F5) 


where ^ n2 (0) is the wave function factor at the y — 0 brane and is given by 
i/v /r +zf cos 2 (7t k n2 )/(2 7T R). The contributions from the folded sector can be extrapo¬ 
lated from the MSSM result by substituting a —> a — 1/2 and Z' v —» Z' Uf . To obtain explicit 
results, we truncate the series and perform a numerical integration over the remaining 
spectrum. The integration bound and the truncation order is fixed such that the numerical 
answer is within < 1% of the analytic answer when Z' v = Z' Vj = 0. 

The left panel of Fig. 14 shows the fractional change of the Higgs soft mass from the top 
sector as a function of (Z' U: Z' Uf ) at a = 1/4. The fractional changes are not very sensitive to 
the choice of a. For comparable Z' v ~ Z' Uf > 0, a < 30% change in the soft mass is observed. 
Corrections to the Higgs mass parameter from wave function factors on the y — 0 brane 
can also be read off of the figure. Noting that at y — 0 the l{ symmetry forces the MSSM 
and folded Z-factors to be equal, it is clear that the corrections from these contributions are 
only 0(10%). 

It is also a logical possibility that the ^-factors are negative. However, in this case the 
transcendental relationship for k n given in Eq. (F4) has a complex solution. We will now 
demonstrate that this leads to the appearance of a zero-norm state in the KK spectrum. 8 As 
a simple example, take a generic bulk hypermultiplet (<F, <F C ) with a non-zero brane kinetic 
term for $ on the y = 0 brane. Assuming the (j) (</> c ) scalar is even (odd) on the y = 0 brane, 


In this respect we disagree with the statement that has been made previously in the literature that Z < 0 
can lead to tachyonic or negative norm states at tree level, e.g. [64]. 
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FIG. 14: Fractional shifts in the Higgs soft mass (left) and term (right) due to brane-localized 
kinetic terms Zy, Z' Uf at a = 1/3. For the purpose of the NDA estimate, we take AR = 5. The 
fractional shifts depend only weakly on a, save at a = 0,1/2 where the shifts are identically zero 
and the fractional changes are ill-defined (also at a = 1/4, where the b^ from the top and /-top 
cancels accidentally at one-loop). 


the presence of the brane kinetic term will induce a kinked wave function: 



kn{y}+ay _ n 

- n k n 
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(F6) 
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Then keeping track of the complex value of k n the normalization is given by 


An, — 


sinh(2 7 i Im(fc n )) 
2 7r Im(/c n ) 


+ 


4 7 tR 


[cos(2 7T Re(fc n )) + cosh(2 7r Im(/c n ))] 


-1/2 


(FT) 


Using the transcendental relation that defines k n , it is straightforward to show that the only 
consistent solution is A n = 0 for the complex k n . Therefore, the spectrum contains a state 
with zero norm at tree level. 

The presence of this state may indicate that negative Z-factors could not be generated 
by well-behaved UV theories. For example, this state should receive a kinetic term at one 
loop; in principle, this can be negative, violating unitarity. In any case, the computation 
underlying Fig. 14 does not obviously generalize to include states with zero norm. We leave 
a full investigation of the impact of negative brane-localized Z-factors to future work. 

We have checked that if one excludes this zero-norm state, NDA-sized negative Z- 
factors can lead to 0(1) corrections to the Higgs mass parameter. Thus the assumption 
of calculability, which is central to the arguments made in the main text, may only hold for 
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positive Z-factors. However, it is unclear whether or not the contribution of the neglected 
zero-norm states would yield a compensating contribution to the Higgs soft mass. 

A similar computation can be carried out for the 6 /t -term. The one-loop diagrams are 
given in Fig. 10, with the result 
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3 y\\i 
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d 4 l 
(2 7r)' 
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n\,ri2 


ni + a + i> n2 (0 ){k n2 + a) 


_ [(ni + a) 2 + l 2 ] [( k n2 + a) 2 + £ 2 ] \ 


+ ( folded sector J 


(F8) 


The folded sector contribution is obtained by substituting a —>• a — 1/2 in that of the MSSM. 
The right panel of Fig. 14 shows the fractional change of the 6^-term coming from the top 
sector at a = 1/4. Again, for comparable Z' v ~ Z' Uf > 0, the change in by is around < 20%. 
The symmetry at y — 0 will result in correction of order 0(10%) similar to the Higgs 
soft-mass. 

This highlights the challenge to calculability posed by brane-localized kinetic terms when 
the 5D cutoff is not parametrically separated from the inverse radius of compactihcation. 
Nevertheless, we conclude that our results are valid across a wide range of parameter space. 
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